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Abstract 

Sleptons can easily be found at future linear e + e~ colliders if kinematically 
accessible. Measurements of their masses and decay distributions would then 
determine MSSM parameters. This paper presents a detailed MC study of 
the production and decay of the lighter scalar tau lepton, f±. We found that 
mfj and 9f (the left-right mixing angle of f) would be measured within an 
error of a few percent. It is also found that tan (3 is determinable in some 
region of the parameter space through simultaneous studies of f\- and e-pair 
production: the polarization measurement of the r leptons from f\ decays 
and the Mi, m^o determination using cr pair production and decay. We also 
point out the possibility to determine e through the measurement of the 
angular distribution of the e/j-pair production. The error on the coupling is 
expected to be comparable to its typical SUSY radiative correction, which is 
proportional to log(mg/mr). The radiative correction affects M\ and tan/? 
determination, necessitating the full ldoop radiative correction to the e~R 
production processes. The implication of these measurements of the MSSM 
parameters on selecting models of the origin of supersymmetry breaking is 
also discussed. 
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I. 1. INTRODUCTION 



The Minimal Supersymmetric Standard Model (MSSM) [TJ is one of the most promising 
extensions of the Standard Model (SM). It predicts the existence of superpartners of SM 
particles (sparticles) below a few TeV to remove the quadratic divergence which appears 
in radiative corrections to the SM Higgs sector. The model is thus free from the so-called 
hierarchy problem inherent in any non-SUSY GUT models. It should also be noted that the 
gauge couplings unify very precisely at high energy in the MSSM, consistent with a SUSY 
SU(5) GUT prediction @. 

Supersymmetry is, however, not an exact symmetry of Nature; instead it should be some- 
how broken to give a mass difference between each particle and its superpartner. Various 
attempts have been made at explaining the existence of soft SUSY breaking Those 
different models of SUSY breaking lead to different relations among the soft breaking mass 
parameters at some high energy scale M SB ; this scale could be as high as M ph or as low 
as (9(10 4 ) GeV, depending on the models. Evolving the mass parameters with the renor- 
malization group equation(RGE) of the model from M S b to the weak scale M weak , one thus 
ends up with different sparticle mass spectra. 

Precise measurements of masses and interactions of superparticles will be one of the most 
important physics targets once they are discovered. If the precision reaches a certain level, 
we will be able to test if a new particle satisfies relations predicted by supersymmetry. It 
will also enable us to measure SUSY breaking mass parameters, and to discriminate between 
models of even higher energy scale responsible for SUSY breaking. 

This is indeed the case at proposed future Linear e + e~ Colliders(LC) operating at ^/s = 
500 GeV , which are designed to provide a luminosity in excess of £ = 30fb^ 1 /year 

f||5|||. It should also be stressed that the background from W boson production to SUSY 
processes can be suppressed drastically thanks to the highly polarized electron beam avail- 
able only at linear e + e~ colliders [|5]J|. 

The production and decay of the lighter chargino, Xi an d scalar leptons, e and jx at an 
LC are studied extensively in previous works f||5|||0,|[ . In particular, it has been shown 
by Monte Carlo (MC) simulations that some relations among soft SUSY breaking parame- 
ters, which are predicted in Minimal Supergravity models (MSUGRA) can be tested very 
stringently It has also been pointed out that one can verify some SUSY relations, 

such as that between the off-diagonal elements of the chargino mass matrix and the mass of 
the W boson, or that between gauge boson-fermion-fermion and gaugino-sfermion-fermion 
couplings J7|. 

In this paper, we discuss production and decay of scalar leptons f and e#, and show how 
various MSSM parameters can be measured from their production only. 

The f is a very interesting object to study, since its mass parameters depend very sensi- 
tively on physics at the GUT scale (Mgut) [0- In SUSY-GUT models the f is in the same 
multiplet with t above M GU t- Therefore the f is expected to have a very large coupling, 
proportional to the top Yukawa coupling, to colour triplet higgs bosons predicted in GUT 
models. Even though all sfermions have equal mass at M p i in MSUGRA models, the large 
Yukawa coupling reduces the f mass at Mgut compared to those of the other scalar leptons, 
which might be regarded as a signature of quark-lepton unification at the GUT scale. This 
observation implies that the stau can be found earlier than the other charged sleptons, which 
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is also phenomenologically interesting. 

In order to obtain the GUT scale mass parameters, one has to evolve the mass parameters 
at M weak toward M GUT . This requires knowledge not only of the f mass matrix at M weak , but 
also of the weak scale r Yukawa coupling Y T = — gm T / (\/2m w cos/5), which is determined 
by the ratio of vacuum expectation values tan/3. Y T may have very big effects on the RG 
running of for tan/? ~ 50: such a large value of tan/3 is expected in a minimal 5*0(10) 
GUT model [lj. 



The measurement of Y T is known to be difficult || , but it has been pointed out |12|JT3 
that the decay distribution of f contains some information on this coupling. 

f production and decay is different from e or ft, because of the non-negligible r Yukawa 
coupling involved in its mass matrix and interactions. Due to this coupling, ?r and fi mix, 
and the mass eigenstates are not necessarily current eigenstates. The same Yukawa coupling 
appears as a non-negligible ttH® coupling, where is a neutral higgsino. This interaction 
is involved in f decay into a neutralino (Xi) an d t, or a chargino(xr) an d v T , since the x's 
are mixtures of higgsinos and gauginos. 

Another feature of f decay that distinguishes it from other slepton decays is that the 
daughter r lepton from the decay f — > further decays in the detector, which enables 
us to measure the average polarization of the r (P T (f — > tx))- The r lepton from fi decay 
is naturally polarized. The polarization P T in the decay f\ — > x1 T depends on Y T . This 
dependence arises because the interaction of gauginos with (s)fermions preserves chirality 
and is proportional to a gauge coupling, while the interaction of higgsinos flips chirality and 
is proportional to Y T . P T from decaying fx reflects the ratio of the chirality flipping and 
conserving interactions and is therefore sensitive to Y T . 

P T also depends on the f left-right mixing angle 9? , and on the neutralino mixing Nj 
which in turn depends on (M 1; M 2 , /i, tan/3). 9? can be determined independently from a 
measurement of the f pair production cross section. On the other hand, information on 
must be obtained elsewhere, for example from e R pair production and decay. Selectron 
pair production involves t-channel exchange of neutralinos. By studying cr pair production 
followed by the decay tR — > ex±, one can thus not only measure the mass of the Xi ( m x?)' 
but also very strongly constrain the gaugino mass parameter Mi. Making use of the mea- 
sured P T {f — > txX) an d assuming a GUT relation between Mi and M 2 , we can in principle 
determine all the parameters of the neutralino mass matrix: M%, [i, and tan/3. One purpose 
of this paper is to reveal the feasibility of the tan (3 measurement at future LC's. 

Another aspect of the ep and f x measurements is also treated in this paper. In the 
high energy limit, cr production involves s-channel exchange of the U(1)y gauge boson 
B and t-channel exchange of its superpartner B. The process turns out to provide clear 
information on the B-cr-cr coupling g B g R£R - Assuming that the e# angular distribution can 
be reconstructed from that of daughter electrons, we find that the sensitivity to the coupling 
9Be R e R wou ld reach (9(1%) (which corresponds to a few % sensitivity to the production cross 
section). The sensitivity is then comparable to the typical radiative correction to the SUSY 
relation g BeR ^ R = V%g', which is proportional to log(m,j/m z ~). This is the first example where 
radiative corrections to couplings involving superpartners might be measured experimentally. 

We also discuss what the P T (f — > TXi) measurement implies in the limit where Xi 
is dominantly gaugino. In this limit, the sensitivity to tan/3 disappears since no ttH\ 
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interaction is involved in the frx\ coupling. However, in this case, we show that sensitivity 
to the chiral nature of frB coupling emerges, offering another test of a supersymmetry 
relation. 

The organization of this paper is as follows. In Sec. 2 we review the physics involved in the 
f\ mass matrix. In Sec. 2.1, we describe the relation between the weak scale parameters m^, 
rrir 2 , and Of and GUT scale f mass matrices in detail. The importance of measuring Of and 
tan (3 is stressed there, since it allows us to check the relation between mg and nif at Mgut- 
Sec. 2. 2 is devoted to describing the procedure to determine tan/? from the measurements of 
P T (f — > 7~x) and of pair production. In Sec. 2. 3, we discuss the energy distribution of f 
decay products, from which P T (f — > rx) and m fl are measured. 

Our MC studies of f pair production and decay are described in detail in Sec. 3, where one 
can find our error estimates on rrif , Of, and P T for J Cdt = lOOfb -1 . Some preliminary studies 



have been given in proceedings reports jl3|, where the effects of the e + e~r + r~ background 
were not properly taken into account. In this paper, we present our final results with an 
optimized set of cuts to remove the background, while minimizing acceptance distortion for 
parameter fitting. These cuts are detailed in Sec. 3.1. The results of the fitting are discussed 
in Sec. 3. 2. 

In Sec. 4.1, we define a function called A% 2 , which allows convenient estimates of errors on 
MSSM parameters that could be obtained through fits of f and e decay distributions. Sec. 
4.2 is devoted to the tan/? determination from a simultaneous fit of f\ and e~R production 
using Ax 2 , demonstrating a unique opportunity to measure tan/3 if it is large. In Sec. 4. 3, 
we go further to determine e (f) coupling to neutralinos. Sec. 5 then summarizes our results 
and concludes this paper. 



2. PHYSICS OF fx 

2.1 Origin of Supersymmetry Breaking and the Mass of f 

tl{r) is the superpartner of tl(r), the third generation lepton. This makes f a unique 
object in the context of SUGRA-GUT models [ID|. 



In minimal supergravity models, SUSY breaking in a hidden sector induces a universal 
soft breaking mass mo, a universal gaugino mass Mo, and a universal trilinear coupling 
Aq through gravitational interactions at the Planck scale M v \. If the soft breaking masses 
remain universal from M p i through Mgut, this boundary condition results in the universality 
at the weak scale of sfermion soft breaking masses within the same representation of the 
SU (3) x SU (2) x U{1)y gauge interactions in the MSSM as long as their Yukawa interactions 
are negligible: 

m\\ weak = (m% UT ) 2 + 0.5(M Gl/T ) 2 (la) 
m\\ weak = (m^) 2 + 0.15(M G ^) 2 . (lb) 



Here is the soft breaking mass of the superpartner of a left(right)-handed lepton, 

rriQ UT and Mq UT are the universal scalar and gaugino masses at Mgut- The model also 
predicts the following relations among gaugino soft breaking mass parameters: 
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M l = ^M^ UT } M 2 = ^M^ UT , M 3 = ^M° UT , (2) 
3a G a G a G 

where M±, M 2 , and M 3 are the masses of U(1)y, SU(2), and SU(3) gauginos (called bino, 
wino, and gluino, respectively). a G is the gauge coupling at Mqut where the aj's unify. 

Eq.(|l|) does not apply for f, due to the possibly large r Yukawa coupling allowed in 
MSSM. The Yukawa interaction of the third generation fermions is described in the MSSM 
by the following superpotential: 

W Y = Y T H\E C L + Y b H 1 D c Q + Y t H 2 U c Q, (3) 

where Hi(H 2 ) is the higgs doublet with hypercharge Y = —1/2(1/2) that gives masses 
to down(up)-type fermions after SU(2) x £7(l)y symmetry breaking: (Hi) = (ttj^Ij 0) 
and (H 2 ) = (0, -^v 2 ). Fermion masses are thus not simply proportional to their Yukawa 

couplings but depend on vi >2 as well: Y T = — gm T j[\f2m w cos/?) where tan/? = v 2 jv\. In a 
simple 5*0(10) model where all the Yukawa couplings are unified at the GUT scale, Y T ~ 1 at 
M G ut and tan/3 is predicted to be around 50 [fTTf] . For such a large tan/3, the contribution 



of the t Yukawa interaction to the RG evolution of mf is non-negligible, vrtf receives a 
negative radiative correction going down from M G ut to M weak , leading to a mass reduction 
compared to m g . Numerical values of the f soft breaking masses at M wea k for a unified 
Yukawa coupling at the GUT scale Y GUT = 1 can be obtained using fbl 



m 2 u = 0.53(M GC7T ) 2 - 0.12m H \ 2 GUT + 0.77m\ GUT (4a) 
m\ = 0.15(M GC/T ) 2 - 0.23m H \ 2 GUT + 0Mm\ 2 GUT . (4b) 

Here m\ G uT is the universal soft breaking mass of sfermions, and tu^gut is the soft breaking 
mass of higgs bosons at M G ut, which may be different from m\ G uT in the S'O(IO) modelJil 
One can subtract Eq.(|l|) from Eq.(||), after setting uihIgut and m\ GUT equal to m GC/T , to 
single out the maximal possible effect of the Yukawa RG running from M GU t to M weak . The 
effect reduces the coefficient of (m? UT ) 2 from 1 to 0.32 for m\ and to 0.65 for m~ . 

u 7 H T J^t 



Yet another source of the reduction of nif is left-right mixing ||15|| . The mass matrix of 
a slepton flavor (II, Ir) can be written as 

M 2 = ( ™Il mill \ = ( m\ + m 2 l + 0.27 D -m z (A + fi tan/3) \ 

\m 2 LR m 2 RR ) \ -mi(Ai + /itan/3) m 2 R + m 2 + 0.23D J ' W 

Here is the higgsino mass parameter, A[ is the coefficient of the soft breaking term propor- 
tional to 1* r -Il-Hi, and D = — m| cos(2/3). The left-right mixing element (m 2 LR ) is negligible 
for the lighter generations. However, for f, if tan/3 ~ 50, the suppression from a factor of 
m T is compensated as long as the diagonal mass parameters are 0(mw)- The mixing is also 
non- negligible if m^ T ,m^ /i. Mixing makes the lighter mass eigenvalue m fl lighter than 
diagonal mass terms. The mass eigenstates and eigenvalues are expressed as 



1 The original formula in Ref. [14] contain D-term contributions, which we have neglected here. 
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cos Of sin Of 
— sin Of cos Of 

,2 i 2 



m LL + ™rr T V (™L - m^) 2 + 4(m| ii ) 1 




(6a) 
(6b) 



tan^f = — - — 2 • (6c) 



m LR 



fx may hence be lighter than e, even in a model with a common soft breaking sfermion mass 
at M weak . 

We learned that determination of tan f3 characterizing the RG running of m? from Mqut 
to M wea k, and of the weak scale f mass matrix parametrized by m^, m? 2 , Of, is necessary to 
extract and at M GU t- The values at M G ut are interesting since they sensitively 
depend on the nature of quark-lepton unification, as has been emphasized recently in Ref. 
|IIJ. The reason is the following: In simple grand unified models such as supersymmetric 
5*0(10) or 577(5) models, the Tr(l) superfield is in the same multiplet as the top quark 
superfield above Mqut- Thus from M p i to Mqut, the t R ( L ) supermultiplet is subject to the 
same Yukawa interaction as the top quark. This reduces ms (and mi T ) at Mqut from that 
at M p i for the SU(5) (5*0(10)) GUT model fpTOfl . The reduction is predicted as a function 



of the top Yukawa coupling Y t , m , M , and A . m R and could be as light as zero at 
Mqut for a large value of A , even if ^ 0. 

Phenomenologically the MSUGRA-GUT suggests that f\ can be the lightest charged 
SUSY particle, thus to be observed first, or might even be the only SUSY particle to be ac- 
cessible at the proposed next generation linear e + e~ colliders. However, we should stress that 
there exist models which predict totally different soft breaking mass parameters 
Dine, Nelson, Nir and Shirman recently constructed a relatively simple model which dy- 
namically breaks SUSY at some intermediate scale [~ 10 6 ~ 7 GeV](DNNS model) H. The 
breaking is then transferred to our sector by a U(l)y gauge interaction, whose scale M m is 
(9(10 4 ) GeV. Its prediction on the gaugino mass parameters turns out to be the same as that 
of the MSUGRA model. This is not the case for the slepton masses, which are predicted to 
be common to (II, v{) and Ir, respectively at M m : 
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*i«§UHs)- (7b) 

Unlike in the SUGRA-GUT model, the slepton masses do not run too much from M m to 
M wea k, as M m is considerably closer to M wea k and there is no strong Yukawa interaction 
involved at these energy scales.@ 

The determination of mi T , rn RT , and mgra at the GUT scale would therfore give us 
a good handle to distinguish the MSUGRA and DNNS models, or if the scale of SUSY 



2 Even if we take a very large Y T at the GUT scale, the weak scale value of Y T is smaller compared 
to that of 1ft and Y t . 
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breaking is below or above the GUT scale; it is not enough to only observe mf < rrie, 
but tan/3 and Of must be measured to determine mi and ms at Mqut- This shows the 
importance of precision studies of production and decay of fx at future LC's. We discuss 
in the next subsections how we can measure these parameters using f\ pair production and 
decay 



2.2 Determination of MSSM Mass Parameters from Production and Decay of 

Sleptons 

Information on Of and tan f3 can be extracted solely from the production and decay of f 



HqJ13|| . In this subsection we sketch our strategy to do this. The determination of will 
be discussed in Sec. 2.3. 

A f decays into a chargino x~ (i=l,2) plus a u T , or a v T plus a W~ or a H~ , or a 
neutralino x°(i=l,...,4) plus a r. Here the neutralinos are some mixtures of the neutral 
components of gauginos and higgsinos (B, W, H® and H® ), and the charginos are some 
mixtures of the charged components. Throughout this paper we assume that the lightest 
SUSY particle (LSP) is the lightest neutralino Xi- Due to R parity conservation in the 
MSSM the LSP is stable and escapes from detection. The decay products of any SUSY 
particle contain at least one Xv 

When both of the pair-produced f\ decay into x?+ i~, the event yields a simple acoplanar 
t + t~ final state. If f\ decays to heavier x° are allowed, the event might contain associated 
jets or leptons. Notice that if the xl, decays into uvxx, the event has the same signature as 
that of f\ — > TXi- If one or both fi'(s) decay into xt + Z/ V the event results in only one or 
zero r lepton + jets or a lepton + missing momentum. 

The fi-to-ino decay branching ratios depend on the scalar tau mixing Of and the pa- 
rameters of the -ino sector (Mi, (M2), /i, tan/3). The measurement of the f\ branching 
ratios might give us extra information on these parameters, but the existence of various 
decay modes also makes the analysis of T\ production very complicated. This point has 
been discussed in previous works |12| in detail, and we will not repeat it here. 

Hereafter we concentrate on the case in which the tr is the second lightest SUSY particle 
and decays exclusively into TXi- Fig.l shows the interactions of the neutral components of 
gauginos and higgsinos with tr and r. The interaction is completely fixed by supersymmetry. 
Namely, the coupling of the tr to B is proportional to the U(l) gauge coupling gi, while the 
coupling to is proportional to Y T . The two interactions have different chirality structure. 
The (super-)gauge interaction is chirality conserving, while the (super-)Yukawa interaction 
flips chirality (In the figure, the arrows next to the f and r lines show flow of chirality). 
Since the polarization of the r lepton P t (tr — ■> txT) measures the ratio of the chirality 
flipping and the conserving interactions, it is sensitive to tan f3. 

As we mentioned already, the gauginos and higgsinos are not mass eigenstates, but 
they mix to form the neutralino mass eigenstates x°(i = 1, ..,4). The ¥r and ?z, also mix. 
Hence the x^Tir couplings depend not only on tan/3 but also on the stau mixing Of and 
the neutralino mixing iV^, where iVy is defined by Xi — N^B + N i2 W + N i3 Hi + N^H^. 
Therefore, the measurement of P T alone can not uniquely determine Y T unless Of and Nij 
are specified. For example, in the limit where the lightest neutralino Xi is a pure bino state 
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(Nil ~^ 1) and m the fi — > tr(l) limit, P T is expressed as 



^ (~ ^ \ 4 sin 2 Of — cos 2 # f .„ . 

P T (n -> 5r = - . ' — - 8a 

v / 4 sin 2 0* + cos 2 6f 



TR - X°it) = V- (, (8b) 

v^iVntan^) + (Y T N 13 ) 2 



(V2Y T N 13 ) -g*(N 12 + N n tajie w )* / x 

^r(r L ^Xl^) = y— ^2 ( 8C ) 

(V2y r Ar 13 ) + (Ni 2 + Nn tan 

respectively. In the gaugino dominant limit, P T does not depend on tan/3 as expected. On 
the other hand, if Ni 3 is non-negligible, P T depends on tan (3, but how it depends differs as 
sin#f varies from to 1. The interactions involving f, Xi xti and the dependence of P T on 
these interactions are listed in appendix A, together with the definitions of neutralino and 
chargino mixing angles. 

Now let us turn to the determination of f and N^. 
• fi Mixing Angle f : 

Since a polarized electron beam will be available at future linear e + e~ colliders, the 
mixing angle Of can be determined by the measurement of the production cross section for 
e + e~ — > fi f{" ||12|| . This can easily be seen by taking the limit mz <C and Pe — +!• 
In this limit, the f production solely proceeds through the exchange of the U(l)y gauge 
boson B. The hypercharge for f L ( R ) is — 1/2(— 1), thus <j(fn) ~ 4cr(fi). Though the cross 
section also depends upon m fl , it can be separately extracted from the energy distribution 
of f decay products or from a threshold scan ]12| . 



• Neutralino Mixing Angles Nj [§,0: 

The neutralino mixing Nij depends on M l5 M 2 ,/i, and tan/3. If we assume the GUT 
relation Mi = haij (3a 2 ) • M 2 , we can determine two out of the three parameters using e^e^ 
pair production as we will discuss below .1 Combining it with the o>f and P T measurements, 
one can then determine all the parameters of the neutralino mass matrix in principle. 

The cr pair production proceeds though the s-channel exchange of gauge bosons and the 
t-channel exchange of neutralinos. We list the amplitudes for the e + e~ — > e^e^ production 
in Appendix B. In the limit Mi, \i 3> mz, \/~s m z and P e = +1, the amplitude reduces to 



iM. — > ifi/g 2 tan 2 Ow sin# 



l-2cos^ / + /3? + 4M 1 2 /s 



(9) 



where and /3f are the polar angle and the velocity of the e^. The first term in the square 
bracket corresponds to the s-channel exchange of gauge bosons, and the second term to the 
t-channel B exchange. One can see that only the interaction with the U(l)y gauge boson 
(B) and the gaugino (B) is relevant in the limit, since e# is an SU(2) singlet. 



3 For any numerical calculation in this paper we assume the GUT relation, though the ratio might 
be determined model independently, using chargino production 
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The differential cross section da(e + e~ — *> e^e^)/ cos# is very sensitive to M\. In Ref. 
H, e^e^j production and their subsequent decays e% — > e ± x± have been studied in detail. It 
was pointed out that the three-momentum of cr can be derived from the momenta of the 
final-state electron pair with a twofold ambiguity, provided that rrig R and m^o are known. 
The cr and Xi masses can, on the other hand, be determined from the energy distribution 
of the electrons with an error of O(lGeV). The study of cr therefore provides two out of 
the three parameters of the neutralino sector. The remaining freedom of tan /3 can then be 
fixed by P T (fi — > r). 

In order to illustrate how the above procedure works, we calculate the cross section 
contours for e + e~ — ► e^e^, and P t (tr — > rx?) 5 fixing m^o at 100 GeV and varying Mi and 
tan/3. Curves in the M-y-fi plane which satisfy the Xi m ass constraint are shown in Fig. 2 
for different values of tan/3. With the mass constraint, one can specify the position in the 
parameter space of the neutralino sector by M x and tan /3, up to a twofold ambiguity of 
positive and negative \i solutions for the large M\ and tan/3 region, or up to a threefold 
ambiguity of two solutions in the negative \i and one solution in the positive fi regions for 
small values of Mi and tan /?.! 

The cjg+g- contours corresponding to positive and negative fi solutions are shown in 

Ft Ft 

Fig. 3a) as the solid and dotted lines, respectively. The difference of the two solutions is 
bigger for smaller tan (3. For tan [3 > 10, the difference becomes negligible. 

The dependence on tan/3 is also mild for <r(e + e~ — ► e^e^), since it only comes in though 
the effect of the gaugino-higgsino mixing, which is suppressed by m 2 z / max (Mi, fi) 2 and 
sm 2 9w, compared to the leading term. The effect is visible for tan/3 < 5 but it essentially 
vanishes for tan/3 > 10. On the other hand, the cross section is very sensitive to Mi 
as expected: it decreases monotonically with increasing Mi, and turns out to be extremely 
small when Mi ~ ^/s, where the t-channel and s-channel diagrams almost cancel each other. 

Fig. 3b) is a contour plot of P T (f R — > x\ T ) i n the Mi-tan/3 plane. As long as Mi is 
not very close to m^, the polarization depends on tan/3 sensitively in the region of the 
parameter space shown in the figure. As one can easily see from Figs. 3a) and 3b), if we 
know Mi precisely from the 6r production cross section, we can extract tan /3 by measuring 
Pt{tr — > Xi T ) unless Mi ~ m^o. Notice when Mi ~ m^o, the lightest neutralino is gaugino 
dominant and there is no significant Yukawa coupling involved in P T as shown in Eq. (|Sa|) . 
Therefore we cannot expect any sensitivity to the tau Yukawa coupling in such a region of 
parameter space. While for Mi 3> m^o, the lightest neutralino is higgsino dominant and Xi 
has significant higgsino component. In such case, some sensitivity to tan/3 is expected for 
a moderate value of tan (3 where the first and the second terms of the numerator of Eq. (|8b[) 
are comparable. 

Notice that in Fig. 3 we did not exclude the region forbidden by the Minimal Supergravity 
Model. In MSUGRA, one has to require the square of any scalar mass parameter be positive 



4 The ambiguities in [i might be removed for tan/3 < 10 by measuring other processes such as 
chargino production and decay. 
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at M p i to prevent the potential from being unbounded from below there. This condition leads 

to the following inequalities at M wea k: m e - fl > ^Jo.87Ml + 0.23D (See Eq.fll])). For instance 
mg R = 200 GeV requires Mi < 215 GeV. If we find Mi > 215 GeV, it will immediately bring 
us to conclude that the SUSY breaking scale M S b is much lower than M GUT , and above 
Msb the theory is different from the MSSM. In the following numerical calculations, we 
will not assume the positivity of the scalar potential at M GU t, since the existence of models 
with Msb "C Mqut is not excluded. The DNNS model is an example of such a model with 
Msb <C Mqut, although their resulting slepton and gaugino masses at the low energy scale 
are consistent with the positive scalar mass requirement at Mg^I 



2.3 Energy Distribution of f Decay Products 

In this subsection, we discuss the measurements of m^, a??, and P T from the decay 
distribution of the r leptons from the f\ decays. As discussed in Sections 2.1 and 2.2, these 



parameters are important to determine both the f mass matrix and tan/5. In Ref. |[L2|| , one 
of us (MMN) has proposed measurements using the pion energy distribution in f — > Xi T ~^ 
x\kv t - In this paper, we discuss measurements using the decay chain f — > Xi r ~^ XiP^t, 
since we find this channel advantageous over r — > x°7ri/ T , as explained below. 

First consider the primary decay f — > rx\- The kinematics is analogous to the e or jx 
cases studied in Ref. ||. The r energy distribution is flat between the endpoints given by 



E* T ±p* T P f 



- C/ max(min) r. ~ ' v / 



where E* and p* are the r energy and momentum in the parent r rest frame: 

m| - m| + ml , 

K = ^ , P* = y/(E* T )*-m* T , (11) 

and [3f = (1 — 4m?/s) 1//2 is the f velocity in the laboratory frame. Knowledge of the two 
endpoint energies allows us to determine mf and m^o, unless (3f is very close to one. 

However, the r decays into Av T , where A = ev e , ir, p, a\, etc., and p^ further decays 
into Ti ±r K° and af to tt ± tt ± tt t or ^n ^. Thus the signature of the f + f~ production is an 
acoplanar two-jet event with low multiplicity. 

In the limit E T 3> m T , the decay products keep the original r direction. However, 
the visible energy is smaller since some of the r's energy is carried away by neutrinos. In 
order to determine mf and m^o , one must reconstruct the original r endpoint energies from 



5 Another important sets of constraints could be obtained from requiring the scalar potential 
neither be unbounded from below(UFB) nor have charge or color breaking(CCB) minima deeper 
than the standard minimum [^6[. One would then find strong constraints on mf L and m,f R at the 
weak scale, depending on fi and mn? H2 . In this paper, we do not consider these constraints, since we 
will not specify mf 2 , w? H , and fi in the later analysis. 
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the energy distribution of the decay products. In Figs. 4a) and b), we show the energy 



distributions of the p and % from a decaying t r(l) with a fixed E T in the limit E T 3> m T |L7 . 
The energy distributions in the cm. frame is obtained by convoluting these distributions 
with the r energy distribution, which we show in Figs. 4c) (for fx — > r^) and 4d) (fi — > r£) 
for a representative set of parameters m fl = 150 GeV, = 100 GeV, and yfs = 500 GeV. 

The it energy distribution depends on P T very strongly. It is harder (softer) for a 7r~ 
from a t^, l ^ (see Figs. 4c) and d)), due to angular momentum conservation. However, a 
substantial correlation between P T and the number of identified events is expected due to 
the inevitable E vis and JPj< cuts to remove the e + e~r + T~ background |]]. As we will see 
in Sec. 3, applying these cuts drastically reduces events in lower energy region E < E^ in , 
where most of the events reside for P T ~ —1. (The maximum and minmum energies of the 
original r lepton -Emax(min) are shown in Figs. 4c) and d).) If E^ in ^/"t*) ^ i s thus hard 
to measure the energy distribution precisely, which results in large errors on E^ in and P T . 
This uncertainty on P T also affects the determination of E^^, as the energy distribution 
near -E^ax depends on P T strongly. Finally, the acceptance depends on P T , giving extra 
uncertainty on the f\ total cross section measurement. 

The p mode is preferable to the n mode in these aspects. The dependence of the energy 
distribution of p mesons on P T is mild, since kinematics forbids low energy p mesons. The 
energy distribution is peaked near E T min for any P T . The P T dependence of the energy distri- 
bution near E T max is also moderate. Because of this pseudo-P T -independence, we can carry 
out the determinations of mf, (m^o), and the cross section without any strong correlation 
to P T . 

Furthermore, the polarization of the p meson depends on P T very strongly, which can be 
seen in the distributions of Pl(t) m Fig. 4c/d) (dashed lines). Namely, a decays mostly to 
a longitudinally polarized p meson (p^) and a t£ decays mostly to a transversally polarized 
p meson (pr)- One can thus determine P T by measuring P p , which in turn can be determined 
from the distribution of the p decay products. A pF decays into 7r ± 7r , and the distribution 
of E n ± in the Pl(t) ~^ tt^tt decay is a very simple function of z c = E n ±/E p , where E p is the 
total energy of the jet to which the 7r =1= belongs, and can be written in the following form 

dT{p T -> 2n)/dz c ~ 2z c (l - z c ) - 2ml/m 2 p , (12a) 
dT{p L -> 27i)/dz c ~ (2z c - l) 2 , (12b) 

where we have ignored terms 0{m 2 p / E 2 p ) but retained O^m^/m 2 ) contributions, and z c is in 
the range (1 - /3 r )/2 < z c < (1 + (3^/2, with V = y/l - Aml/m? p . 

By fitting the z c distribution together with the E- ]et distribution, one can determine both 
P T and m,f x . The error from the small P T dependence of the E- )et distribution is reduced by 
the simultaneous use of the z r distribution. 
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3. MONTE CARLO SIMULATIONS 



3.1 Event Selection 



In this section, we investigate the feasibility of the f studies outlined above at future 
linear e + e~ colliders. 

As discussed in the previous section, measuring the E p distribution of the cascade decay 
fi^r^p and the z c (= E n ±/E p ) distribution of the subsequent p decay p — > ^n , one 
can determine m^, P T , and Of- However, in order to measure these parameters, one has to 
introduce cuts to control backgrounds. We also have to reconstruct p, a%, or 7T from f\ de- 
cays with minimum mis-ID probability among these channels. The reconstruction efficiency 
heavily depends on detector performance, necessitating Monte Carlo (MC) simulations with 
realistic detector and machine parameters. In this subsection we dicuss the dominant back- 
ground from e + e~ — > e + e~r + r~ and present our cuts and detector set up to reduce it. 
We also define our cuts to identify p and a x and MC-examine the contamination due to 
misidentifications. The result of the fit to MC data after the cuts will be presented in Sec 
3.2. 

In the following we study a sample case of fi rf pair production followed by exclusive 
fi — ► TXi decays. We will not treat the other decay processes: — > T^Xii > 2) or v T % i 
where the expected event signatures are much more complicated. The helicity amplitudes 
for r^ff production and their subsequent decays into rx? are calculated using the HELAS 
library [TE|]. The final state r leptons are generated using the BASES/SPRING package 
|l9fl , and are decayed with TAUOLA version 2.3 [ 20[ . The effects of initial state radiation, 



beam energy spread, and beamstrahlung are also taken into account |ZT] . 

The end-product stable particles (7r ± , 7, e, ) are then processed through a detec- 
tor simulator, and are identified, if possible, as (7r =l= ,7...) candidates. In this paper, we 
assumed the JLC1 detector parameters, except for the forward electron veto system. The 
model detector is equipped with a central drift chamber (CDC), electromagnetic and hadron 
calorimeters(EMC, HDC), and muon drift chambers, whose parameters can be found in Ref. 
f| . The used detector simulator is the same as the one used in the previous studies |J . 

The event signatures for the f pair production are acoplanar two jets or one jet + one 
lepton. The former mode is cleaner since the latter mode suffers from WW, ei/W, and 
eeWW backgrounds. We will, therfore, concentrate on the former mode. We used the 
following basic cuts to select such an acoplanar two-jet event: 

1. There exist two and only two jets for some y cnt > 2.5 x where y cnt is imposed on 
the reduced jet invariant mass: EiE 2 ■ (1 — cos6> 12 )/(-E v is) 2 > Vcut- 

2. Both of the two jets must clear the polar angle cut | cos 6^1 < 0.8. 

3. The net charge of each jet must be unity and opposite in sign to that of the other. 

4. The acoplanarity of the two jets has to be large enough: 6^,,^ > 30°. 

5. These two jets have to have invariant masses consistent with r hypothesis: mj et i, 
"^jet2 < 3 GeV. 
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6. The missing transverse momentum (Pt) has to exceed 15 GeV. 



7. There has to be no electron or position above 6 e > 50 mrad from the beam axis. 

In addition we need cuts to identify r decay products as p or a% in order to analyze each 
decay mode separately: 

• p cuts 

A jet with two 7- and a ^-candidates is identified as a p^ — > 7r ± 7r° candidate if 
m 27 < 0.25 GeV and m jct < 0.95 GeV. 

If there is only one 7 candidate in the jet (7 + 7r ± ) , we require m- ]et < 0.95 GeV, 
assuming a possible cluster overlapping in the calorimeter. 

• a\ cuts 

A jet is identified as a\ if it contains three charged 7r's only, or four or three 7's + one 
7r , or two 7's + one 7r =!= with m2 7 > 0.25 GeV. 

Cuts 1) through 4) are similar to the one used in the previous studies ||. These cuts 
together with cut 5) were designed to reduce the background from gauge boson productions 
(WW, evW, eeWW) to less than 1 fb for P e = 0.95 at yfs = 500 GeV. 

Cuts 2), 4), 6) and 7) are to remove the e + e~ — > e + e~r + r~ background, where the two 
photons radiated off from the initial-state e + and e~ collide to produce a r pair, while the 
e + and e~ escape into the beam pipe. Most of proposed LC detectors have a relatively big 
acceptance hole in the forward region. For example the JLC1 model detector can only veto 
e s above 6 > 150 mrad, allowing the two-photon background with a f*? up to as high 
as 75 GeV kinematically, and up to 37.5 GeV typically, since it is quite rare that the two 
initial-state particles give the maximum possible transverse kick in the same direction to the 
t + t~ system. 

In the previous Monte Carlo studies of the backgrounds for sfermion productions, the 
e + e~l + l~ backgrounds were eliminated by cuts on acoplanarity angle and Pt @- The same 
applies to the eerr background in the f production studies. One might worry about the 
r decay giving extra JPj< to the eerr events, however, the overall reduction of the energy 
of the r decay products compensates this effect. The JP? > 35 GeV cut together with an 
electron veto angle of 9 v e eto = 150 mrad, the cuts on polar angle (cut 2), and acoplanarity 
angle (cut 4) turns out to remove most of the background events. We have generated 110K 
e + e~r + r~ events with 9 e < 150 mrad, E T1 + E T2 >15 GeV, | cos# r | < 0.9, and > 10°, 
using a code developed by Kuroda f22[.i The corresponding production cross section is 1.10 



pb, and therefore the generated events are about / Cdt = lOOfb equivalent. 12 events 
survived cuts l)-5), 7) and J> T > 35 GeV: a(eerr)| cut = 0.12fb±0.035fb. 

Introducing such a high missing Pt cut might introduce an extra correlation between 
the acceptance (or the measured value of the f pair production cross section) and the 



We have generated the e + e _ r + r~ events in the phase space suffiently larger than the one defined 
by the cuts 2) and 4), because the reconstructed jet axes do not in general coincide with the original 
r directions. 
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polarization of r lepton from f\ decay, because jet energies become softer for a smaller 
P T (f — > tx) as discussed already in Sec. 2. 3. Our MC simulation for m fl = 150 GeV, 
m^o = 100 GeV, and ^/s = 500 GeV, shows that 17.3%(P T =1) /17.1%(P T = -1) of the 
generated signal events are identified as fx with no JPt but P v ; s > 10 GeV cutQ, while 
only 12.6%(P T =+1)/9.8%(P T = -1) of them identified for JP T >35 GeV.i Though the 
correlation between P r and the acceptance is smaller than that of the 7r mode, and P r can 
also be constrained from the z c distribution of the p decays, the reduction of the acceptance 
by up to a factor of 2, and its strong P T dependence might be worrisome. (See Fig. 5 for 
the /Pr distribution of the signal events corresponding to 10 4 generated f pairs and the 
eerr background for / Cdt = lOOfb -1 .) It should also be noted that events with smaller jet 
energies are less likely to be accepted. This might complicate the simultaneous measurement 
of m^o and rrif, as the measurement of the energy distribution near P^m becomes more 
difficult (see Fig. 6). 

In this paper, we therefore assume a forward coverage down to 50 mrad. The e + e _ r + r~ 
production cross section for E T1+T2 > 15 GeV, | cos# T | < 0.9, 9 acop > 10°, and 9 e < 50 mrad 
is 0.719 pb. Out of 70K generated e + e~r + r~ events , which correspond to / Cdt = lOOfb -1 , 
only 19 events remained as background after applying cuts 1-7). The overall detection 
efficiencies for the signal events after the same cuts are 17.0% and 16.0% for P T = 1 and 
P T = —1, respectively. Jet energy distributions of p-identified events for different P-r cuts 
are shown in Fig. 6 for = 150 GeV, m^o = 100 GeV, and P T = — lB One can see that 
the Pt > 15 GeV cut has no significant effect on the signal events as expected, while for 
the Py > 35 GeV cut, the acceptance diminishes drastically for Pj et < 50 GeV, making the 
determination of E^ in difficult. 

In order to realize the 50 mrad veto angle, we need to place additional veto counters in the 
beam background mask, which might produce extra beam backgrounds. On the other hand, 
having a tighter forward veto can significantly reduce the SM background in the low P T 
region, which will help us extend our discovery reach to SUSY particles with a mass which 
is very close to that of the LSP. Further studies are necessary to optimize the parameters 
for the extra forward electron veto. Another possibility to reduce the Pt cut value is of 
course to go down close to the f pair production threshold. If the f production is accessible 
at -^5=350 GeV, fx > 25 GeV must be enough to eliminate the eerr background. 

Now we are going to discuss the p and a\ cuts. These cuts are chosen to minimize 



7 The acceptance is smaller than that of the other sleptons since we had to require both of the 
r's decay hadronically. The background to the search mode where one t decays into e/p is larger 
but expected to be manageable. We also had to apply a tighter jet polar angle cut to reduce the 
e + e~r + r~ background. 

8 The P T dependence comes mostly from f\f\ — > irp events, where 7r's tend to have low energy for 
P T = —1. Those events are less likely to be accepted due to the Pt cut. 

9 The jet energy distributions are slightly softer than that of Fig.4d), as the MC simulation includes 
beam effects and initial state radiation. These effects will also be included in the fits of Sec. 3. 2. 
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contaminations of p to a\ or vice versa due to mis-reconstruction of photons. As described 
earlier, the p and cti decays involve 7r°'s, which in turn decay into 27's. For a high energy 7r°, 
however, the two photons are occasionally misidentified as a single photon due to the cluster 
overlapping in the calorimeter, and therefore the a\ sometimes has the same signature as 
p. Fig. 7a) shows the jet invariant mass distributions of the events consisting of a tt~ and 
one or two photon candidates coming from f + f~ — > t + t~ — ► 7i + p~ and 7r + a^. The solid 
histogram is for p~ decays from 50K f pairs forced to decay into 7i + p~, and the bars are 
the data of 7190 f pairs forced to decay into 7r + a|[. The ir + p~ data are scaled so that its 
relative normalization to n + ai is correct. Due to the mis-reconstruction the events have a 
considerably smaller jet invariant mass distribution compared to that of g^'s decaying into 
2tt~tt + (dotted histogram). The p cut on the invariant mass mj et < 0.95 GeV only removes 
half of the eti contamination. 

In Figs. 7b) and c), we plot the jet energy and z c distributions of the p candidates that 
satisfy the p cuts in the same 7r + p~ and 7r + a^f samples. Solid histograms are those from 
the 7r + p _ sample and bars are from the 7r + a^ sample. The number of identified p events 
from 7r + p~ I V + ai samples is 4522/400 for P T = +1. The contamination is larger for a higher 
Ejet and a lower z c . This is because high energy 7r°'s from a\ decays have less chance to be 
identified as 2 photons, thereby sneaking into the p signals. The same MC simulation told 
us that very few a\ decays could be reconstructed with iV 7 > 3 if Ej Ct > 50 GeV. 

The contamination affects and P T fit to E )et and z c distribution. In principle the a\ 
contamination to the p sample must be corrected for before the data are fitted to obtain 
rrir or P T . However, due to rather low expected statistics (~ 1400 p candidates expected to 
survive after the cuts for 10 4 f pairs), we did not attempt making such corrections at all. 
We will see in the next subsection that input parameters of MC and the corresponding best 
fit values of the m f and P T fits to O(10 4 ) fi pair events are consistent with each other. 



3.2 Fit to MC data 

In this subsection we present the results of our fits to the selected MC data for represen- 
tative sets of parameters, with m fl = 150 GeV, m^o = 100 GeV, and y/s = 500 GeV, and 
backgrounds for / Cdt = 100 fb" 1 and P e = 0.95. 

In sections 2.1 and 2.2, we have discussed the importance of measuring and Of to 
determine the weak scale f mass matrix, and P T to determine tan (3. These parameters can 
be measured by looking at p candidates from fi cascade decay (Sec. 2. 3): m fl is measured 
through the energy distribution, sin^ through the production cross section a flfl) and Pi- 
through the z c distribution. 

In the following, we fit the E- ]et and z c distributions of the p candidates selected from 
the signal MC data to numerical functions calculated by convoluting the r — > Pl(t) decay 
spectra with the r energy distributions. The fit parameters are m^, m^o, P r , and the 
number of produced fi pairs N flfl . The results of the fit to 10 4 and 5 x 10 3 fifi pairs will 
be shown in this section. Notice that the production cross section of tr(tl) with rrtf = 150 
GeV is about 0.09 pb (0.02 pb) for P e = 1; therefore, the generated 10 4 fx events roughly 
correspond to / Cdt = 100 (400)fb~ 1 , respectively (see Fig. 1 of Ref. HI~2] ). The fit will be 



15 



extended to include the measurement of e# production and decay in Sec. 4, to obtain the 
error on tan f3. 

We first describe our calculation of the theoretical distributions for the fit. As we men- 
tioned earlier, the Ej et and/or z c distributions are numerically calculated by convoluting the 
r — » Pl(t) and the r energy distributions.^ In the calculations of the theoretical distribu- 
tions, we took into account the effect of beamstrahlung and initial state radiation on the 
fi energy distribution as before pi|. We also found that the acoplanarity angle cut (cut 



4) has a significant effect on the energy distribution. Since the acoplanarity angle cut is 
very complicated to implement in the numerical calculation of the energy distribution, we 
approximated the effect by imposing an acolinearity cut of 30 degrees in the CM frame of the 
f\ pairs instead. The resultant jet energy distribution that was calculated this way roughly 
reproduced the shape of the energy distribution of a statistically larger MC event sample 
(10 5 f\ pairs). The overall normalization has been determined by comparing it with the MC 
simulation, and corrected for a P T -dependent acceptance factor. Agreement between the z c 
distribution of the MC data and that of the numerical calculation was poor for z c ~ or 
z c ~ 1 due to the acceptance effects; we therefore fit the z c distribution only over the range 
0.08 <z c < 0.92. 

In the previous subsection we have seen that the Pt > 15 GeV cut is necessary to reduce 
the eerr background for an electron veto of 50 mrad. Since the signal events are hardly 
affected by the cut, and since it is hard to implement the Pt cut in the numerical calculation 
of the fitting curve, we decided to ignore the Pt cut for both the MC events and the fitting 
curve, and neglected the eerr background altogether, though the dominant backgrounds 
from WW, eeWW, ZZ, and vvZ corresponding to J Cdt = 100 fb^ 1 have been included in 
the fits in this section.0 The production cross sections for the backgrounds before and after 
the selection cuts and the number of remaining p events are listed in Table 1. 

We first separately perform a fit of and m^o to the Ej et distribution, and that of P T 
to the z c distribution. Figs. 8 a) and b) are the results of our mass fit to the Ej et distribution 
of 10 4 fifi events decaying into trXi exclusively. 1476 events were identified as p and used 
for the jet energy fit. In this fit, we kept P T = +1 and set the normalization of the curve 
so that the total number of events agreed with that of the MC data. Fig.8a) plots the jet 
energy distribution of the MC events together with the best fit curve obtained by minimizing 
the log-likelihood function (= x 2 ) with and m^o as free parameters. The contours for 
Ax 2 = 1 and 4 in m^o-m^ plane are shown in Fig 8b). The MC events were generated 
for m fl = 150 GeV and m^o = 100 GeV, while the best fit values are m fl = 147 GeV and 
m^o = 97 GeV for this MC sample. The values are consistent with the inputs; Am^o=2.8 
GeV and Am^ = 3.9 GeV can thus be expected as 1-cx errors on these quantities. 

The errors on the two parameters might be reduced further. Notice that we have only 



10 We included the effect of the finite p width for the jet energy distribution as in Ref. [17] but did 
not take it into account for the z c distribution. 

11 £^ v i s > 10 GeV is implicit for all the MC event generation in this subsection. This condition is 
not included in the fitting curves as its effects were found negligible. 
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used the events identified as f\ — > r — > p here. The other modes into tt, a±, or other leptons 
can also be used to increase the statistics. One may combine the information from other 
sparticle decays once they are observed. The previous analysis of jl decays showed that 
Am^o = 1 GeV can be achieved typically, which would reduce Am^ down to 1.5 GeV. 

On the other hand, P T can be determined by fitting the it fraction of the parent p 
energy (z c = E^/E^t). In this fit, we fixed m^o and to their input values and the 
normalization of the fitting curve to the total number of the MC sample. Figs. 9a) and b) 
show the z c distribution for the selected p candidates together with the best fit curve. Here 
we used the data in the region 0.08 < z c < 0.92 and £j et > 20 GeV since the low energy 
region is insensitive to P T . 924(885) MC events were used for the fit where fi's decayed 
into 7k(t£,)'s. The best fit values and their errors were obtained to be 0.995 ± 0.082 and 
—0.991 ± 0.08 for tr and tl, respectively. 

In order to justify such separate fits of E- ]et and z c distributions with some of the fitting 
parameters fixed by hand, we must make sure that there is no strong correlation among m fl , 
m^o, P T , and cr^^. For this purpose, we calculated the errors on the masses and P T by fitting 
a two-dimensional distribution in (Ej et ,z c ), varying m^o, m^, P T , and the total number of 
produced stau pair events N TT (which corresponds to if the integrated luminosity and 
acceptance is known). 1224 events in the interval 0.08 < z c < 0.92 were used for the fit with 
P T = 1. 

The resultant errors obtained from this fit agreed very well with the previous estimates. 
The best fit value for the masses are = 146.3 GeV and m^o = 95.4 GeV. The shape 
of the x 2 contour projected onto the m^o-m^ plane looked quite similar to that of Fig.8. 
The estimated errors are Am^ = 4.07 GeV and Am^o = 2.99 GeV, being consistent with 
the previous estimates, taking into account the difference in the number of events used for 
each fit. The error on P T was also calculated allowing m^o, m^, and N TT to move freely 
in minimizing \ 2 - The best fit value of P T is P T = 0.89 ± 0.07, again consistent with the 
previous estimate. These results support our assumption of small correlation between the 
energy distribution and P T for the r — > p decay mode. 

Finally we move on to the determination of the f mixing angle Of. As discussed already, 
our strategy is to use the measurement of the production cross section together with that of 
rrif . We generated 5000 f\ pairs decaying into tx\ with P T = 0.6788. The SM backgrounds 
for / Cdt = 100 fb _1 were also included. The p signal sample corresponds to sinflf = 0.7526, 
and m? x = 150 GeV with a bino dominant LSP. 

We used events in the region 0.08 < z c < 0.92, which, after selection, reduced to 628 
events. Since the acceptance differs by about 12% between P T = 1 and P T = —1 due to the 
z c cut,0 we minimized \ 2 in the N flfl - plane, varying P T , m fl) and rn^o. No significant 
correlation was found between and m^o and the estimated errors are Am^ = 6.6 GeV 



The data in z c ~ or 1 can of course be used once detector performance is understood and 
included in the numerical calculation of fitting curves. The dependence of the acceptance to P T 
described here is purely artificial, unlike that caused by Pt cuts. 
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and Acr^ = 2.2fb.lI 3 J The mass error is consistent with the previous estimate if one takes 
into account the difference in the numbers of produced events and the acceptances. The 
error on the cross section is consistent with the error simply estimated by the statistics of 
the accepted events. 

Fig. 10 plots contours of constant minimized y 2 surfaces projected onto the m^-sin^f 
plane. We found Asinftf = 0.049. One can see that the correlation with makes the 
error large. It is possible to reduce the f\ mass error by using the m^o obtained from other 
measurements as it was discussed earlier. Am^o = 1 GeV would reduce Am ?1 to 2.21 GeV, 
in which case the error on sin# f is less than 0.03. However, the error cannot be less than 
0.014, which is limited by the observed number of events.0 

The fx decay process studied above is quite complicated compared to that of e or p. 
Nevertheless, in the above discussions, we have found the measurements of masses m^o, 
P T , and N flfl could be done without any significant correlations each other. We have also 
learned the mass errors using (9(600) accepted p events are consistent with those of (9(1400) 
events if the latter is statistically scaled. The error on the cross section is consistent with 
the error estimated by the statistics of the accepted events. This allows us to estimate the 
errors on m^, P T , Nff, and Of reliably by simple statistical scaling of each error in a wide 
region of parameters space. This fact is used when we combine the f\ measurements with 
the e measurement in the Sec. 4. 

In this subsection we assumed that a tight forward electron veto is possible and applied a 
small Pt cut as was discussed in the previous subsection. If this Pt cut has to be increased 
to 35 GeV, the result of this subsection must change. The mass measurement is based on 
the measurement of -Emax(min) extracted from the energy distribution of p candidates. As 
can be seen from Fig. 7, events near E^ ay[ are not affected by the cut, therefore AE^^ will 
not change, either. E^ ax is sensitive to m fl — m^o, therefore Am fl obtained with m^o from 
other slepton measurements will not change significantly. AP T will not be affected too much, 
either, since the z c distribution is not sensitive to P T for E p < 20 GeV, which is the region 
most affected by the Pt > 35 GeV cut. Finally, some extra dependence of the acceptance 
on P T should be introduced by the large /Pt cut. This might increase the error on the 
production cross section, since the acceptance moves by 20% with P T varying from 1 to — 1. 
However, the dependence can be tamed by measuring P T from the z c distribution. 



4. COMBINED ANALYSIS 

4.1 Ax| anc * Ax~ functions 

In this section we are going to extract tan/3 by combining the measurements of 
P T (f — > 7~x), sin^f, and the knowledge of the neutralino mass matrix obtained from the 
measurements of production and decay. Some MC simulation of e production had already 

13 We assumed / Cdt = 100 fb -1 . 

14 Including other decay mode to the analysis would improve statistice, and reduce A6 T . 
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been carried out for a specific set of parameters 0, and we have just finished a corresponding 
MC analysis for f\ in the previous section. It is now straightforward to perform combined 
fits to determine the MSSM parameters for representative points in the parameter space. 
Nevertheless it is quite time consuming to do it exactly as in a real experiment. 

Therefore in this paper, we define Ax 2 -like function Ay 2 to estimate the sensitivity of 
LC experiments to these SUSY parameters. As we mentioned already, we take a sample 
case where both fi and e# are produced at a future LC. The e# pairs are selected requiring 
acoplanar e + e~ paris, while the f\ pairs are taken from an acoplanar two-jet sample. Since 
these samples are statistically independent of each other, we define the e and f\ parts of the 
Ax 2 functions A%~ and A%| separately in this subsection. A%~ and A%? are functions of 
two sets of input parameters: (mg R , Mi(M 2 ), //, tan /3) and (m~ , M((M 2 )', //, tan /3') for Ax~, 
and (mf 1 ,9f,M 1 (M 2 ),ii, tan/3) and (m' fl , d' f , M[(MJ), //, tan /3') for A*?. The A^ 2 function 
are defined in such a way that Ay 2 = when the two sets are equal and the projection of the 
hypersurface of Ay 2 = 1(,4, 9..) to one of the parameters (m~ , M[....) fixing (rrie R , M 1 ...) 
roughly agrees with the 1(,2,3..)— a error of that parameter. In this subsection, we first 
define A%§,-n in detail, and then discuss the error on tan (3 in the next subsection. Readers 
who are interested only in the results can skip this subsection. 

The polar angle((9) distribution and the end point energies of electrons from e# decays 
can be measured at future LC experiments as discussed already in Sec. 2.2. Therefore we 
define A\ 2 ~ by using the two sets of quantities as follows: 



Axl(me R , M 1 (M 2 ),fJ,, tan/3; m' gR , M[(M 2 )\ //, tan/3') 

K -nif 



"bin („,l _ ,, \2 




(13) 



where rii and v! i are the expected numbers of events in i-th bin between —1 + 2(i — V)jnun < 
cos 9 < —1 + 2i/nbin calculated for the first and the second sets of input parameters: (mg R , 
Mi(M 2 ), /i, tan/?) and ijn'^ R) M((M 2 )', //, tan/3'), respectively. For later use we calculated 
rii assuming / Cdt = 20 (or 100) fb _1 , ^/s = 500 GeV, a 27 % acceptance, and n bin = 25, 
making use of formulas for the pair production cross section listed in Appendix B. The 
acceptance is chosen to be a factor of 0.6 smaller compared to the value obtained by the MC 
simulation ||. This is because our expression for the selectron production cross section does 
not include any effects of initial state radiation, beam energy spread, and beamstrahlung. 

E^ x and are the upper and lower end points of the energy distribution of electrons 
for the first (the second) set of parameters. AE^ ax and AE^ in are defined by 



A£ma x( min) = J X E Hn , (14) 



ay 



where iV av = n tota iE bin / (£^ ax - E^ in ) with E bin = 4 GeV being a kind of bin width, and 

^total = J2i=l n i- 

Ax? is chosen to reproduce the actual A\ 2 of the MC data fitted with the second set 
of parameters (m' gR , M[(M 2 )', //, tan/?') when n i; iV av ^> 1 and n tota i ^> iV av . The reason 
is the following: If r^'s are replaced by actual data, the first term of Eq.([13|) is \ 2 °f the 
data fitted with the parameters (m' gR , M((M 2 )', tan/3') based on Gaussian distributions. 
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The difference between the data and n«'s divided by Hi must be small in the limit of large 
statistics if (mg R , Mi(M 2 ), /x, tan/?) is the parameter set that nature has taken, therefore 
the projection of the hyper-surface that satisfies A%| = 1 to one of the fitting parameters 
roughly indicates the size of ±l-cx deviation of the parameter from the best fit point. In 
this sense, we can call the first set of parameters as input parameters and the second set as 
fitting parameters. 

The second and the third terms are intended to represent the sensitivity of the electron 
energy distribution measurement to determine -E^ax an d E^ in . As has been mentioned 
already, mg R and m^o can be determined from the endpoints of the energy distribution of 
electrons. In actual experiments, electron energies are measured by some detector with a 
finite energy resolution. The JLC1 model detector, for example, has <Je/E = lh%/\fE © 
1% ||, demanding us to take into account the finite size of E hin ~ 2 GeV. Moreover the 
energy distribution can also be smeared and distorted due to finite beam energy spread, 
beamstrahlung, initial state radiation, and possible dependence of acceptance on electron 
energies. In the previous study the shape of the energy distribution and its dependence on 
mg and m^o were obtained from the MC study itself, and the mass errors were estimated by 
actually fitting the energy distribution. 

For simplicity we assume here that the energy distribution is flat between -E^ax and E^ in , 
while conservatively taking Eun— 4 GeV. If the average number (iV av ) of events in a single 
bin is large enough so that the fluctuation is negligible compared to N av , the central value 
of -E max (min) is obtained as 

-^max(min) = -^cnd ± ^ ^ E hin , (15) 

where E^ nd is the central energy of the upper (lower) edge bin, N hin is the number of events 
in it, and iV av is the average number of events in some intermediate bin. Based on the 
statistical error on iVbin estimated assuming a Gaussian distribution, the error on -Emax(min) 
is A£^ lin(max) /£ , bin = 1/V-Wbin- 

However, when the actual -Emax(min) * s vei T c l° se to a bin boundary, the fluctuation of 
iVbi n becomes non-Gaussian. iVb in can even exceed iV av or becomes zero, making nonsense 
of Eq . (p~5|) and the error estimate. Therefore in Eq.(|l4|), we replaced N^n by N av /2 , which 
corresponds to the choice of binning where -Emax(min) * s approximately at the center of the 
edge energy bin. In such a case, the fluctuation of the edge bin becomes Gaussian-like as 
long as iV av is large enough, thereby justifying our estimation. 

Several comments are in order. One might think that the measurement of the end point 
energies does not fit any \ 2 analysis implying a Gaussian distribution, if the energy resolution 
is too good and the expected numbers of the events in the edge bins are too small. In such 
a case the probability distribution for -E^ax or E^ in is expected to be asymmetric, because 
if an event is observed in some energy bin between E\ and E 2 , selectron and neutralino 
masses which give E^ ax < E\ or E 2 < E^ in are strongly disfavoured. However, the number 
of events in a single bin is expected to be large enough in the "precision measurement" phase 
of the LC , and hence our treatment assuming a Gaussian distribution can be justified: we 
have checked if our treatment of AE 1 ^^^ roughly reproduces the previous results on ft 
production and decay 0, and found that the A% 2 contours by using the last two terms of 
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Eq.(|13D but replacing -Emax(min) to -^max(min) i n the m /r m x? pl ane agree very well with the 
previous results for the same number of accepted events. 

When we calculated the 9 distribution or the electron energy distribution we occasionally 
found rii less than 15, or 3Ai<^, min N > E hin /2. In such a case, we merged the cos# bins or 
enlarged E^. Our treatment underestimates the sensitivity compared to any log-likelihood 
analysis based on a Poisson distribution and therefore is conservative. 

The Ax 2 analysis mimics the true x 2 fit to the 9 and E e distributions, though it neglects 
the correlation between the 9 and E e distributions through the total number of events: the 
fluctuations of the events in E e bins have correlations with the fluctuations of rtj's, since the 
events must add up to an equal number in both distributions. This correlation disappears, 
however, in the limit where the number of events in the edge E e bins is negligible compared 
to the total number of events, thereby justifying our method. 

Finally, in the definition of Ax|, we assumed that the e production angles are recon- 
structed precisely. This is not true since there is always a wrong solution of 9 for each 
event. The wrong solutions must first be included in the 9 distribution, which must then be 
subtracted statistically, bringing more uncertainty into our analysis. We also assumed that 
all the selected selectrons contribute to the determination of the production cross section, 
angular distribution, and masses, and will not distort the measurement due to decays 
into heavier neutralinos. 

For fx pair production and their cascade decay t\ — > r — ► p, we have already discussed 
that it is important to measure the total cross section, P T from the decay distribution of the 
p decay products, and E^ oin and -E^ax from the energy distribution of p candidates. Fits to 
10 4 fi pairs have been done in Sec. 3.2, taking into account backgrounds corresponding to 
/ Cdt = 100 fb _1 , and the errors on m fl , m^o, and N flfl have been obtained. In the following 
we define the f\ part of our Ax 2 -hke function Ax? so that it reproduces the results in Sec. 3. 2: 

= —ff— + { AE^ ) + {-^—) + {-AIT) ' <16) 

where is defined to be the sum of constant background (iVbg) and the total number of 
signal p events (iV t J al ) for which both of f's decay directly into Xi T an d the r's then decay 
hadronically. We took N^ g = 100. iV to ( al was estimated using an integrated luminosity of 100 
fb _1 , the tree level cross section without any beam effects, and the acceptance obtained in 
the previous simulation with no fx cut. The branching ratio to rx? was calculated by the 



formula in Ref . |12j . Notice that in the region where the lightest neutralino is higgsino-like, 
the lighter chargino Xi an d the lightest and the second lightest neutralinos x? an d X2 are 
almost mass-degenerate. Then the decay modes f\ x°r and f% — > x v t generally open 
up, which would yield rather complicated final states with associated jets. As we have not 
studied the sensitivities and backgrounds to these modes, we will not include them in the 
study below, instead we will simply take the number of events where both f\& decay into 
Xi to estimate N^ tal - 

The first term of Eq . (|16|) is intended to show the statistical significance of the total fi 
pair production cross section. On the other hand the second and the third terms express 
sensitivity to -E^ax an d E^- m . We again calculate AE^ DSSL , J]An \ using a rather simple set of 
formulas: 
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AE^ X = AEZ ta = l.Sy/2/N^, (17) 

where iV av = EunN to tai/ (-^max ~ -^min)- The effect of the smearing of the energy distribution 
by the cascade decay of f\ is taken into account by the overall factors in the right hand 
side of Eq.(|T?D. The factors are chosen so as to reproduce Am^ and Am^o in the previous 
subsection. The larger factor for AE^ ax compared to AE^ in may be understood as the effect 
of the higher reduction of the events near E^^ due to the r decays( see Figs. 4c), and d)). 
Finally, AP T is estimated statistically scaling the error on P T in Sec. 3. 2: 



AP T = 0.07 x (1400/^/1400 + N hg ) / (n^/Vn) . (18) 



In Eq.(16), we assumed that there is no large correlation among the measurements of 
-Emax(min)' an d ^r- We also imply that errors on these parameters can be estimated by 
the statistical scaling of the results in Sec. 3. 2. These features have been checked explicitly 
by the MC analysis in the same subsection. 



4.2 Determination of tan j3 from slepton production 

We have already pointed out in Sec. 2. 2 that the simultaneous measurements of fi and 
€r productions would determine tan (3. In order to estimate its statistical error expected at 
a future LC, we have defined A%?/ f ) functions in the previous subsection. A%| is a function 
of two sets of MSSM parameters (m^, Mi....) and {ml- ,M[....). In the limit of infinite 
statistics, the projection of the hypersurface of A% 2 = 1(,4,9..) to one of the parameters 
(m^, M[....) fixing (m e - H , Mi...) agrees with the 1(, 2, 3..)-<r error of that parameter obtained 
by using the cos 9 distribution and the end point energies of electrons from e.R decays of real 
data. The definition of Ax% is similar, but here the data used are the number of signal p 
events from fx — > Xi T followed by r — > u T p, the end point energies of the r's the fi decays, 
and the average r polarization. 

We will start our discussion with the determination of the parameters of the neutralino 
mass matrix from e# production alone. m e - fl and m^o are determined essentially by through 
the energy distribution of the electrons from decays. On the other hand, Mi is mainly 
constrained by the e# production cross section. The dependence of the total cross section on 
Mi and tan/3 for fixed mg R and rn^o, assuming the GUT relation between Mi and M 2 , has 
been shown in Fig. 3a). Notice that constraining mg R is very important for the determination 
of Mi, as the production cross section depends not only on Mi but also on rrie R . 

Fig. 11 shows the error on Mi estimated with the A%| function, where the input pa- 
rameters were chosen such that m^o = 100 GeV, P e = +1, mg R = 200 GeV, p > 0, and 
/ Cdt = 20 fb _1 . Both positive and negative errors on Mi are shown in %. The errors were 
calculated by finding a minimum value of Ax 2 for a fixed M[ = M 1 + AM 1 varying the other 
fitting parameters (m'g, p', tan/3'), with the MINUET program. The values of AMi which 
give Axmin = 1 are plotted as 1-cr lines in the figure. This corresponds to projecting the 
Axj = 1 hypersurface to the Mi axis. Mi will be determined within an error of 5 ~ 7%, 
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typically in the region of parameter space shown in the figure .113 

We found that the errors are asymmetric for positive and negative fluctuation of M\. 
This is because the condition A% 2 = 1 forces the total cross section to be close to its input 
value. The cross section increases with decreasing tan /3 for tan /3 < 5 and decreases with 
increasing Mi for fixed m^o, as can be seen in Fig. 3. For tan/3 = 1.5, the error on Mi is 
therefore larger in the negative fluctuation (Fig. 11): Increasing tan /3 can compensate the 
increase of Off due to the reduction of Mi here. For tan /3 — 15, a similar argument shows 
that the error is larger in the positive fluctuation. 

In the figure, we have also shown the error on Mi for tan f3 = 15, but restricting tan f3' > 
10 for searching minimum value of A%g. In this case, the la fluctuation is symmetric and 
smaller. This suggests that even a rough estimate of tan/3 can greatly help us restrict Mi. 
As we will see below, such improvement is indeed possible in some region of parameter space 
if fi production is observed. 

Now we turn to the determination of tan/?. As we have discussed already in Sec. 2. 2, 
tan /3 can be extracted from the polarization of r leptons produced in fi decays, if we know 
Of and the neutralino mixing angles iVy. 

Of is determined from the measurement of and fi production cross section. (The 
result from our full MC analysis is in Fig. 10 in Sec. 3.2). The sensitivities to the production 
cross section, m^, and P T are taken into account in the definition of the Ax? as the first 
term, the second and the third terms, and the fourth term of the right-hand side of Eq. (|l^) , 
respectively. 

On the other hand some information on neutralino mixing can be obtained from selec- 
tron production, e.g. by minimizing Ax~- The Neutralino mixing depends on only 3 (4) 
parameters Mi(M 2 ), /i, tan/3, where A%? strongly constrains two of these three parameters, 
Mi and m^o, as described earlier in this subsection. 

Notice that the constraint on m^o from Ax? is stronger than that from Ax?, because of 
the larger statistics fot the e production and the smearing effect of the end point energies 
for -E'max(min) > which have been taken into account in the definition of AE^ ax and AE^ in 
in Eq. (|TTD . The small Am^o from the cr production helps determine tan/3 better, for the 
following reason: Am fl correlates with Am^o, and AOf with Am fl , as can be seen in Fig. 8b) 
and Fig. 10. Therefore AOf is smaller for smaller Am^o, which reduces the uncertainty coming 
from f mixing in determining tan/? from P T . 

Fig. 12 plots Ax 2 = Axl + Axf = 1 contours projected onto the Mi-tan/? plane for input 
values of (Mi(GeV), tan/3)= (219.0, 15), (149.3,15), (124.5,15), (124.5,25), and (124.5,5). The 
other input parameters are common for all representative points,mg = 200 GeV, = 150 
GeV, m^o = 100 GeV, and sin 6^ = 1 (fi = tr) and we took / jCdt = 100 fb _1 . 

As we have discussed previously, P T depends sensitively on tan /3 if m^o <C Mi or Xi is 
higgsino-like, because in this case the \fr coupling involves the r Yukawa coupling. Thus 



15 We found that our error on Mi is smaller than previously quoted [||. The previous estimate 
did not use the GUT relation between Mi and Mi and the constraint on m^o which would have 
been obtained from the cr production was not exploited in the Mi and M 2 fit [p3|| . Our result is 
therefore consistent with the previous one. 



23 



the error bar is expected to be smaller for a larger M\. However, the number of accepted 
events becomes small when m^o <C M\ and f\ also decays into X2 or Xii therefore the error 
on tan/5 for M l = 219 GeV is larger than that of M x = 149.3 GeV. When tan/? = 15, the 
lower(and upper) bounds of tan/3 for the A% 2 = 1 contours are 13.85 (18.5), 13.28 (16.37), 
8.94 (18.74) for Mi = 219, 149.3, and 124.5 GeV, respectively. 

Finally, in our definition of A% 2 we did not include any luminosity error. In Ref. || it 
has been argued that the luminosity can be measured with an error of (9(1%). On the other 
hand, using a typical e production cross section 0.1 pb~ 0.2 pb (See. Fig 2), and assuming 
a constant acceptance of 27%, we can see that the errors on M x are estimated based on 
(9(1000) accepted e events for the luminosity of 20 fb -1 (Fig.ll), and (9(5000) events for 
100 fb -1 (Fig. 12). This corresponds to an error on the cross section of about 3% and 1.4%, 
respectively. The latter is already comparable to the luminosity error. Hence a further 
increase of statistics would not improve the actual error on the MSSM parameters unless 
the error on the luminosity is also reduced further. 

The estimated errors on tan (3 are rather impressive, compared to those from the other 
experimental methods. In Fig. 13, we plot the error on tan/5 which can be obtained from the 
lighter chargino (xf) production and the co-production of Xi an d X2 ■ Here use has been 
made of the direction of pair produced charginos xf with a 10% acceptance from its decay 
products, assuming f Cdt = 100 fb -1 for both P e = 1 and 0.lll The errors on m-± and m-± 
were both assumed to be 2% as long as they are accessible kinematicallyHI One can see that 
the upper bound practically disappears when tan/? exceeds 5.0 

4.3 Checking Supersymmetry Relation 

So far we have been assuming that new particles found at an LC are superpartners 
of leptons. In other words, we have implicitly been using supersymmetric interactions of 
sfermions with the neutralinos of the MSSM without any attempt at checking the nature of 
the interactions. Instead, we merely used the data to determine the free parameters of the 
MSSM, such as M x and tan (3. In this subsection we are going to discuss the possibility to 
probe the gaugino-sfermion-fermion interaction (more specifically, the B-e^ -cr coupling), 
and some aspects of the -B-fj-T coupling. 



The direction of a produced chargino can be solved for with a two- fold ambiguity when the 
chargino decays into Wxi ||- The forward backward asymmetry for the final state W can also be 
used even if the xt decays into W*Xi [00- 

Am-i was found to be around 5% for 50 fb -1 of data ||. A threshold scan for the xfXi P arr 
production might determine the mass better. 

18 It has been claimed that a very precise measurement of tan [3 is possible when tan (3 ~ 4 |7j if 
the chargino mass errors are negligibly small. Some additional error on tan f3 has been introduced 
here assuming a finite error on m ±. In Fig. 13, we have also taken larger value of rap compared 
to 0, where sensitivity to tan/3 is smaller. 
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We start our discussion with Sr. The production proceeds though the s-channel exchange 
of gauge bosons and t-channel exchange of neutralinos, whose cross section is shown in 
Appendix B. The t-channel exchange is dominated by bino-like neutralino exchange, which 
led us to the simple dependence of the cross section on the gaugino mass M 1 as has been 
shown in Fig. 3a). 

The tree level coupling of the B-cr-cr vertex has a simple relation to the B-e-e coupling 
in the MSSM: 



9§. 



eReR 



V2gtan6 w = V2g'. (19) 



This relation is imposed by supersymmetry. Thus the measurement of g B - e will allow us to 
prove that e and B are indeed superpartners of tR and B. 
For this test we modify the relation of Eq. (|i~9"D as 



9 B . 



£ReR 



V2g'Y B . (20) 



and estimate the sensitivity to Y B by introducing a new A% 2 function for the selectron pair 
production which depends on Y B though gRe R e- m the limit of mz <C Mi and /i, we obtain 
an approximate formula for the matrix element M.: 



M. cx sin 9 



1 



4Y, 



2cos9[3 f + [3 2 f + 4Mf/s_ ' 



(21) 



It is apparent from Eq . (|2l"D that one can constrain both Y B and Mi by measuring the 
differential cross section: da(e + e~ — > GrSr) / dcosO. 

Figure 14 is a Ax? contour plot projected on the Mi-Y B plane for a representative point 
in the parameter space of the MSSM: m- eR = 200 GeV, /x = 300 GeV, M x = 99.57 GeV and 
tan j3 = 2. One finds a good sensitivity to the coupling Y B of (9(1%) in this case. The reason 
why we got upper and lower bounds on Mi and Y B is as follows: when we increase M[ from 
Mi to Mi + AMi, the total cross section decreases. The corresponding increase of Ax? can 
be compensated by increasing Yr correspondingly. However, due to the constraint on m^o 
which comes from the electron energy distribution from decaying e^'s, the optimized value 
of // is smaller for a larger M[ (see Fig. 2 for the relation between Mi and fi). Hence the 
second lightest neutralino mass is lighter for a larger M[, and it also tends to have a larger 
mixing with B. At the M[ upper bound of A% 2 = 1 , the polar angle distribution changes 
its shape so that it is less forwardly peaked. On the other hand, the lower bound on M[ is 
determined by m^o. Namely, for given Mi, there is a upper bound on m^o which one can 
obtain by varying /i and tan/3 maximally. It becomes thus harder to reproduce the input 
m^o as we decrease M[. 

Some deviation of Y B from its tree level value is expected if we take into account the effect 
of radiative corrections in the framework of the MSSM. If there is a large difference between 
several soft SUSY breaking mass parameters, such a correction occurs very naturally. For 
example, if rrig 3> mj and m^, the effective theory below Q < rriq is not supersymmetric, and 
couplings related by supersymmetry start to run differently according to the RG equations of 
the effective theory. In particular, both squarks and quarks decouple from the wave function 
renormalizations of gauginos in the low energy effective theory, while only squarks decouple 
from that of gauge bosons, from which Y B ^ 1 may originate. 
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The RG equations below the squark decoupling follows 
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da' 55 /2 

-a , (22a) 



dhgQ 12vr 

da B~ee 11 / / 99M 

dhiQ = T, a a ^ (22b) 

where we neglected terms proportional to (Y B — 1) on the right-hand side of the equations. 
We find from Eq.(^) 

A^ee/^ee = 0.007**, (23) 

where t~i = log 10 (m$/m f ). 

It is rather striking that the error on the couping is of about the same order as that 
of the radiative correction proportional to log(m,j/m z ~) if the squark mass is much heavier 
than the slepton mass. This, on one hand, requires knowledge of rrig and a full 1-loop 
calculation of the process to remove the uncertainty in Yg from the determination of Mi, 
notice that the error on M\ increases by a factor of two, if we let Yb move freely. This 
also implies a larger error on tan/3, as the errors on Mi and tan/3 are correlated strongly 
when the lightest neutralino is gaugino-dominant. On the other hand, we can turn this 
argument around. Then emerges the possibility to constrain the squark mass scale from the 
measurement of Yg or other couplings even if the energy of future colliders is not enough 
for the squark production. A full calculation of 1-loop radiative corrections to this process 
is eagerly anticipated. 

A similar radiative correction to gy /iieL turns out to be of the order of 2%x \og 1Q (mg/mj), 
but the sensitivity to this coupling has been argued to be rather poor f7j: about — 15%+30% 
for a representative parameter choice. This estimate is based on the study of gaugino- 
dominant chargino production and decay, using the forward-backward asymmetry of the 
decay products and the total production cross section. The chargino production proceeds 
through t-channel exchange of v, and the v is assumed kinematically inaccessible in the 
study. Its mass is determined by comparing the production cross sections for polarized 
and unpolarized electron beams, and the decay forward-backward asymmetry, but it has a 
very large uncertainty. Furthermore, the uncertainty in the branching ratios introduces a 
systematic error to the measurement of the total cross section. Notice that the production 
cross section of gaugino-dominant charginos is very small for a right-handed polarized beam, 
hence the chargino study heavily relies on the use of the unpolarized beam, where large 
backgrounds limit the decay modes to study. The estimated error on the cross section is 
about 5% for / Cdt = 100 fb -1 . These uncertainties limited the gwp e studyB 

The coupling g B e R e * s considerably easier to measure than g-fype- nas a sizable produc- 
tion cross section for the right-handed electron beam, and the uncertainty in the production 
cross section is expected to be very small if the e decays exclusively into e\\. For the rep- 
resentative parameters of Fig. 14, the t-channel exchange of neutralinos is dominated by the 
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During this work, we learned of similar work on measuring decoupling effects by H.-C. Cheng, 



J.L. Feng, and N. Polonsky \26\. We thank Jonathan Feng for bring this work to our attention. 
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lightest one, Xi > an d its mass is well constrained from the cr decay data. It is therefore not 
surprizing that the coupling is measured very well here. 

Deviations from tree level MSSM predictions can also appear in other couplings involving 
sleptons or neutralinos, such as g^. Unfortunately, the measurement of tan/3 has a large 
error, thus radiative corrections may not be relevant in this case. The gaugino mass matrix 
also gets radiative corrections to its tree level value |27j. If one assumes a unified gaugino 
mass at the GUT scale, one may in principle extract the squark mass scale from the gaugino 
mass relation at the weak scale. Unfortunately, the measurements of gaugino masses are 
limited by ambiguities in the neutralino and chargino mixing angles 0^,0] • The gluino mass, 



though involving no mixing, is hard to measure precisely at hadron colliders, too. [pq 

Now we turn our attention to f\ production. In the previous subsection, we found that 
P T (fi — > TXi) becomes independent of tan/3, if Xi is gaugino- dominant, or in other words 
m^o ~ Mi. In such a situation, simultaneous measurements of Of (using the total fi-pair 

production cross section) and P T constrain the nature of the B-f-r coupling instead of 
constraining tan/3. Given the fact that Xi is almost a pure gaugino (which can be checked 
with scalar electron production), the measurement of the total f\ pair production cross 
section essentially fixes the polarization P T through Eq.(|8a|). Any deviation of P T from it 
indicates that something unexpected is happening. 

In Fig. 15 we show A% 2 = 1 contours by taking the mixing angle parameter (9?) in the 
Xjt coupling free from that in the Zff coupling (Of). We can say that Of — Of measures 
the chirality flipping part of the B(W)-f-T interaction which is zero in the MSSM. Due 
to the dependence of P T on tan/? through a small but finite higgsino component in the 
neutralino mass eigenstate (xi), the sensitivity of the frx? coupling to the Of is marginal 
unless \l > Mi.0 

Nevertheless the figure can be regarded as an example of a "no lose theorem" of precision 
measurements of supersymmetry processes. Depending on the position in the parameter 
space of the model, we occasionally lose sensitivity to some parameter, as we have seen 
for the tan/3 determination using slepton production. We can, however, turn this into an 
advantage: the process becomes independent of the ambiguity caused by the parameter, and 
we can test its super symmetric nature. In the current case, we can check the chiral nature 
of the gaugino, thanks to the insensitivity of the process to tan f3. 

5. CONCLUSION 

In this paper we presented an extensive study of the production and decay of the lighter 
scalar tau lepton f\ at a future LC and discussed physics that could be extracted from them. 
Studying f production is important because it may be lighter than the other sleptons, and 
could thus be found earlier. The light f\ case is also theoretically well motivated in the 
MSUGRA-GUT model, and is not excluded, at least, in other models, as long as there is a 
large f L -f R mixing. 



20 The contours in Fig. 15 depend sensitively on the region of tan j3 searched. We took 1 < tan j3 < 
50 here to obtain Ay; 2 = 1, assuming Yukawa coupling is not too large at GUT scale. 
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We discussed that the f mass matrix at Mgut might provide a clue to distinguishing 
SUGRA-GUT from DNNS models. In order to obtain the GUT scale mass matrix, one 
must know the mass matrix at the weak scale and the tau Yukawa coupling (Y T ) which is 
characterized by tan (3. 

The mass matrix can be determined if one knows the f masses and mixing angle Of- 
The feasibility of determining those parameters at an LC has been studied for the lighter 
mass eigenstate fi, assuming the JLC1 model detector. For a representative parameter set: 
m,f x = 150 GeV and m^o = 100 GeV, we found that these masses can be measured to Am fl = 
4.1 GeV and Am^o = 3 GeV for 10 4 fx pairs produced with background corresponding 
to J Cdt = 100 fb _1 , assuming that fi's decay exclusively into Xi T - F° r the same mass 
parameters and luminosity conditions, the expected statistical error on the mixing angle 
turned out to be Asin^f = 0.045, when sinftf = 0.75. 

The polarization (P r ) of r leptons from fi decay is sensitive to tan/3 because of its 
dependence on the r Yukawa coupling. The expected statistical error on the polarization 
was estimated to be ±0.07 for 10 4 fi pairs and background corresponding to / Cdt = 100 
fb _1 . Using the information of the neutralino mass matrix obtained from the simultaneous 
studies of the cr production and decay, tan/3 might be determined. The error on tan/3 
varies drastically with Mi and m^o, as shown in Fig. 12 for some representative points in the 
parameter space of the MSSM. 

Notice that tan/3 is one of the most important parameters that determine the Higgs 
sector of the MSSM. At the same time, it is known to be difficult to measure especially if it 
is large. If tan (3 > 10, fi decays give us a unique opportunity to determine tan /3. 

We have also discussed a possibility to test the supersymmetry relations among couplings 
involving superpartners. By studying the polar angle distribution of e# production, one can 
measure not only Mi, m^o but also the gaugino-selectron-electron coupling g B g Re - A ^ 
allowing g BiRe to move freely from the tree level prediction of supersymmetry gives Ag B - Re 
~ (9(1% ~ 2%). This is comparable to typical radiative corrections to the same coupling 
~ 0.7% x log(mq/mr). This suggests that the LC might allow us to start probing radiative 
corrections to couplings involving SUSY particles. 

Implications of the MSUGRA model at LEPII and LHC have been discussed and studied 
in many papers. Unfortunately, prospects to determine the soft SUSY breaking mass param- 
eters are not so bright there: as for LEPII, its available luminosity is too low for the slepton 
study and one has to fight the enormous background coming from W + W~ production for 
chargino study ||24j| , while at LHC, one suffers from the high QCD background even though 
strongly interacting superparticles will be copiously produced. Therefore those studies in 
the framework of the MSUGRA model are focused mostly on the discovery potential of 
the machine in question. However, it is becoming more and more recognized that we can 
certainly go beyond that if a next generation linear e + e~ collider is actually built. Namely, 
the experiments at the LC will make it possible to measure the parameters of the MSSM 
once a superparticle is discovered, which will then enable us to check the predictions of the 
models of SUSY breaking. 
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APPENDIX A. PARAMETERS OF THE MSSM 



In this subsection, we are going to summarize the interactions of f that are relevant 
for the analysis in this paper. The interactions are fixed by supersymmetry and gauge 
symmetry, as well as by the mass parameters of the model. 

The f-r-ino interactions relevant to fj decay are expressed by the following Lagrangian 

(24) 



C= Y, f l f{P L af j + P R a^)x]+ Yl nUrPRbijx} + h.c, 



i=l,2 j=l,..,4 



i=l,2 .7=1,2 



where 
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V2m wCOS ^ ^ a Lj = -^[N j2 + N n t&nO w ], 



a% = -^NjitanOw, a Rj 



gm 7 



bij = -gUji, 



N. 



\f2m w cos (3 

bR 3 - ~7= a U &- 



J'3j 
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\/2mw cos P 

Here the real orthogonal matrix and unitary matrices C/y and are the diago- 
nalization matrices of the neutralino mass matrix M.^ and chargino mass matrix M.q as 
follows: 



D 



I ' M ( \ ' 1 = Mq, NMnN- 1 = M N , 
where the mass matrices are written in the following form: 

M N {B,W 3 ,H U H 2 ) = 

/ Mi —m z sinOw cos/5 m z sin#jy sin/3 \ 

M 2 m z cosOw cos/3 —m z cosOw sin/? 

—m z sm.0w cos/3 m z cosOw cos/3 —Li 

\ m z sm.O w sin/3 — m z cosOw sin/3 — \i / 



(26) 
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Tfl\v y/2 COS/3 Li 



(27a) 
(27b) 



29 



fj, is a supersymmetric higgsino mass parameter, while M\ and M 2 are the soft breaking mass 
parameters of Bino and Wino introduced previously. The mass eigenstate Xi an d current 
eigenstates B,W,Hi,H 2 are related by 



N i2 W + N i3 H, + N iA H 2 . 



(28) 



Unlike the notation of Haber and Kane 



], we take N to be real so that m^o can be 



either positive or negative. Its sign must be kept to understand the equations in Ref |fL2| . We 
take Im-col < \m^o\ < \m^o\ < \m^o\ and < m~- < m^-. We assume the mass relation of 

I Ai — I X'2 1 — 1 A3 I — I A4 — X\ — X2 

MSUGRA Mi = (5/3) - tan 2 9y^M 2 for numerical calculations in order to reduce the number 
of parameters. 

Eq.pl) leads to an expression for P T : 



X?r) 



— (a 



L \2 

11) 



(afi) 2 + (a^) 2 



(29) 



APPENDIX B e R PRODUCTION 



cr production proceeds through t-channel exchange of neutralinos and s-channel ex- 
change of gauge bosons. The tree level couplings of the e^eXi vertices may be read off from 
Eqs.(2~3) and fl2"5|) by setting sin#f = 1 and replacing r — > e. We obtain the formula for the 
e^-pair production cross section as follows: 



da (h e ,h e ) = ^^-j:\M(h e ,he)\ 2 



dcos9 



2s 16tt2 



iM(h e ,h e ) = -iXiC^ sin 6s/3 f 



(30a) 
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e 2 (1±(- 

+ 7 + ^2- 
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\ - 1 A 2 jR 
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(30b) 



where h e (h e ) = ±1/2 represents the helicity of the initial-state electron (positron), Aj = 
h e — h e , 9 and /3/ is the e R production angle and velocity, and t = — | (1 — 2cos9j3f + 0f). 
The couplings Ah e and Ajr(l) are given by 



A 1/2 = sin 2 9 W , A_i/ 2 
Aj R = —\Z2gtan9wNji. 



1 

2 



-- + sin 2 9 W 



(31) 
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FIGURES 



FIG. 1. Interactions of neutral components of gauginos and higgsinos with tr and tl or tr. 
FIG. 2. m^o = 100 GeV contours in the M\-u plane: solid and dotted lines correspond to 

Al 

tan/3 = 1.5 and 30, respectively. 

FIG. 3. a):<7-+-- contours with m~ e = 200 GeV, Js = 500 GeV, and P e = 1 in the Mi-tan/3 

e .R e .R 

plane. At each point of the figure, a is chosen so that m^o = 100 GeV. Solid lines correspond to 

Al 

a p > solution and the dashed lines to p < 0. b): P t {tr — ► tx?) contours in the Mi-tan/3 plane 
with the same neutralino mass constraint. Only the contours of p > solutions are shown. 

FIG. 4. a) and b): Energy distributions of the p and 7r from a) tr decay and b) tl decays with 
a fixed E T (E T » r) as functions of z = E^i p \/E T . c) and d): Energy distributions of the p and 
7r from a cascade decay of a f for = 150 GeV, m^o = 100 GeV and y/s = 500 GeV. The f 
decays exclusively into t R( ^ in c)(d)). 

FIG. 5. Pt distributions of events passing cuts l)-5), for yfs = 500 GeV, = 150 GeV, 
and m^o = 100 GeV. The solid and dotted lines with higher Pj> tails are for 10 4 f pairs decaying 

Al 

exclusively into t£ and r^, respectively. The /V distribution of the eerr background, also shown 
in the figure, corresponds to J Cdt = 100 fb _1 . 

FIG. 6. -Ejet distributions of 10 5 f\ pairs decaying exclusively to for different total Pt cuts: 
no p T cut (dotted), p T > 15 GeV (solid), p T > 35 GeV(dashed). p T > 35 GeV is the optimal 
cut for 9 v e eto = 150 mrad, while P T > 15 GeV for veto = 50 mrad. 

FIG. 7. a) invariant mass distributions of jets consisting of 7r~ + one or two 7's for f + f~ 
— ► r + r~XiXi> followed by r + r _ — > n + p~v T v T (solid) and 7r + a^^ r ^ T (bars). The latter corresponds 
to misidentified a{ — > 7r _ 7r 7r — ► 7r _ 47. We assumed the parameters of the JLC1 model detector. 
The invariant mass distribution for the jets from a\ 3-prong decays (a 1 — > 27r~7r + ) is also shown 
as the dotted line. Due to the photon mis-measurements the jet invariant mass distribution of a\ 
1-prong decays sits below the one for 3 prong decays, b) The jet energy and c) z c distributions for 
the events that satisfy p cuts. The solid line is the distribution for f + f~ — > ir + p~ , while bars are 
of f + f~ — > 7r + a^ as before. The contamination is larger for small z c or higher £j et . 

FIG. 8. Results from the mass fit to 10 4 fifi pair events decaying into trXi exclusively, where 
the SM background corresponding to / Cdt = 100 fb _1 has been included in the fit: a): the jet 
energy distribution for the p events selected from data MC events (bars) and the best fit histogram. 
In the fit we kept P T = +1 and normalized the histogram so that the total number of events agreed 
with that of the MC data. The average SM background is also shown in the figure, b) Contours 
for Ax 2 = 1 and 4 in the m^o and plane. 
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FIG. 9. z c distributions for the p candidates selected from 10 4 f\ pairs decaying exclusively 
into &):tji and b):r^ , together with the best fit histogram. The sample with E- jCt > 20 GeV and 
0.08 < z c < 0.92 are used for the fit. The best fit values of P T and their errors were obtained to 
be 0.995 ± 0.082 and -0.991 ± 0.008 for t r and t l , respectively, for fixed m fl = 150 GeV and 
m,o = 100 GeV. 

Xi 

FIG. 10. Ax 2 contours in the m^-sinOf plane, resulting from the fit to 5000 f\ pair generated 
for JCdt = 100 ftr 1 , m fl = 150 GeV, m^o = 100 GeV, and P T = 0.6788. The MC sample 
corresponds to a f\ with sin6> T = 0.7526 decaying exclusively into a bino-like lightest neutralino. 

FIG. 11. 1-cr error band on M\ estimated using the Ax 2 function. Solid and dotted lines plot 
values of M[ — M\ in %, where Axflmin = 1 for the M{ varying other fitting parameters 
fi' , tan/?'). Input values are chosen so that m^o = 100 GeV and /i is positive for each M\. Solid 
lines are for tan/3 = 1.5 as input where the parameter region with 1 < tan/3' < 100, and 50 
GeV< n' < 10 4 GeV was searched to obtain AM^ Dashed lines for tan/3 = 15, but only the 
parameter region with tan/3 > 10 are searched to obtain these 1-cr limits on M\. 

FIG. 12. Ax 2 = A%! +AX~ = 1 contours projected onto the Mi-tan (3 plane. Projections of the 
contours on M\ or tan /3 corresponds to 1-cr errors of the parameter. Input values are chosen to be 
mfj = 150 GeV, sin#f = 1 (j\ = tr), m^o = 100 GeV( p > 0), and m e - s = 200 GeV. Input values 
of n are explicitly shown in the figure for individual sample points. P t (tr — ► rx?) = 0.8,0.4,0, 
and —0.4 contours with m^o = 100 GeV and fi > are also plotted in the figure as dashed lines. 

FIG. 13. 1(2)-'<t' errors on tan (3 from chargino production as functions of input tan/3. We used 
chargino distributions for P e = +1 and with / Cdt = 100 fb _1 , and assuming that both chargino 
masses are known to 2% accuracy. Upper bound practically disappears when tan/3 exceeds 5. 
Input values are M 2 = 210GeV, ^ = -195 GeV and mt, 2 = 500 GeV. 

FIG. 14. Axl = 1 contour in the Mi-Yg(= fl£ g / V2 g') plane. The definition of Ax| has 
been modified to allow g^- e to deviate from y/2 g' . Input values are rrie R = 200 GeV, jjl = 300 
GeV, Mi=99.57 GeV, and tan/3 = 2. The error on the coupling is of about the same order as that 
of the radiative correction proportional to log(mq/m z ~) when mq/mj ~ 10. 

FIG. 15. Ax 2 = 1 contours when the mixing angle Of in the Xi^ r is allowed to move freely 
from that in the Ztt coupling (Of). Of — Of parametrizes the chirality flipping part of B(W)- t-t 
interaction, which is zero in the MSSM. The solid line corresponds to fi = —600 GeV and the 
dashed line to /i = —200 GeV. The other parameters are fixed to Mi = 104.5 GeV, and tan /3 = 20. 
The sensitivity is moderate for /i S> Mi. 
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TABLES 



process 


CJp e=+ o.95(fb) 


c C ut (fb) 


number of p candidates/ (100 fb 1 ) 


WW -► T + T~ 


6.23 


0.16 


14.9 


eeWW r+r" 


2.16 


0.21 


17.9 


ZZ — > t + t~^P 


4.88 


0.59 


51.0 


z/PZ — > r + r~ 


0.46 


0.07 


6.3 



TABLE I. Dominant background cross sections at P e = +0.95 to the process e + e~ — ► f\f\ 
followed by f\ — > x§r. Background cross sections after requiring the cuts described in the text, 
and the average number of p background events for / Cdt = 100fb -1 are also shown in the table. 
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